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Optimizing plant control values of a power plant 

5 



DESCRIPTION 

10 

Field of the Invention 

The invention relates to the field of process control, in particular to a 
15 method and computer program product for optimizing power plant control 
values and a power plant optimizing system as described in the preamble of 
claim 1 , 7 and 8, respectively. 

20 Background of the Invention 

Methods and systems for optimizing power plant operation are shown in US 
5,347,466 and US 5,886,895. In both cases, plant operation, which is con- 
trolled by control values generated by a control system, is optimized in order 

25 to minimize fuel costs while achieving predetermined required output values 
for produced power and process steam. This is done by, at a given time, de- 
termining future values of control values and simulating a behavior of the 
plant up to a given future time. From the simulation, fuel costs are deter- 
mined and a objective function comprising the fuel costs and costs for buy- 

30 ing power from another source is computed. In order to determine optimal 
control values, the control values are varied and the simulation is repeated 
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until a minimum of the objective function is found. However, such an opti- 
mization is constrained by the fact that predetermined operating limits of 
the plant may not be exceeded. This has the effect that the optimization 
does not necessarily give a financially optimal result for the plant owner. 

5 

Description of the Invention 

It is an object of the invention to create a method and computer program 
10 product for optimizing power plant control values and a power plant opti- 
mizing system of the type mentioned initially that overcomes the deficien- 
cies mentioned above. 

s 

These objects are achieved by a method and computer program product for 
is optimizing power plant control values and a power plant optimizing system 
according to the claims 1 , 7 and 8. 

c 

In the inventive method and computer program product for optimizing j 
power plant control values and a power plant optimizing system, a rate of 
20 ageing of plant components is determined when simulating a future behavior } 
of the plant, and an objective function to be minimized by the optimization 
comprises said rate of ageing. 

In this way, the invention incorporates long term financial effects of compo- 
25 nent ageing into short .term decisions on plant operation. For example, un- 
der certain conditions and when energy prices are high, the optimization 
shall indicate that it is financially advantageous to operate the plant outside 
normal conditions, i.e. that higher revenue from generated power more than 
compensates for higher maintenance or replacement costs due to increased 
30 component degradation. 
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In a preferred variant of the invention, plant output values representing plant 
production are determined in the simulation and a difference between said 
plant output values and demand values is included in the objective function, 

5 Further preferred embodiments are evident from the dependent patent 
claims. 



Brief Description of the Drawings 

10 

The subject matter of the invention will be explained in more detail in the 
following text with reference to preferred exemplary embodiments which are 
illustrated in the attached drawings, in which: 

Figure 1 shows a block diagram corresponding to the invention; and 
15 Figure 2 shows a block diagram corresponding to a preferred embodi- 
ment of the invention. 
The reference symbols used in the drawings, and their meanings, are listed 
in summary form in the list of reference symbols. In principle, identical parts 
are provided with the same reference symbols in the figures. 

20 

Detailed Description of Preferred Embodiments 



Figure 1 shows a block diagram corresponding to the invention. An optimi- 
25 zation module 1 is provided with demand values d{r) . The optimization 
module 1 generates power plant operating parameters or control values u(r) 
that are input to a simulation module 2. The simulation module 2 comprises 
a simulation model of an actual power plant 3 and determines values for a 
power plant output or production p(r), revenue values q(r) from selling the 
30 produced output, costs c(r) associated with a given power plant output, and 
values of a rate of ageing e(r) of plant components. The values determined 
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by the simulation module 2 are input to the optimization module 1 . 

The values d{T\u{r),p{T\q{r),c{t\e(T) mentioned above are, in general, vec- 
tors that change over time, i.e. each vector represents a trajectory in a mul- 
5 tidimensional space. For example, 3(t) represents an changing amount of 
electric power and of process steam to be generated over the next hour or 
day. 

The control values «(r) represent plant input or operation parameters that 
10 are used to control the plant, such as valve and guide vane positions, or set 
points for local controllers of temperature, fuel flow or steam flow. From a 
trajectory of control values w(r) that is given for a time interval from a pres- 
ent time / to a future time t + T, where T is called "prediction horizon" or 
"optimization horizon", the simulation model determines a trajectory of the 
15 state of the power plant, i.e. the values of a plurality of variables that repre- 
sent the state of the plant as it changes with time. From the state trajectory, 
trajectories p(T),q(r),c(T),e(T) for the values output by the simulation module 
2 are determined for said time interval. In a preferred embodiment of the in- 
vention, the simulation model is a dynamic model of the power plant and its 
20 components. 

The production p(r) represents the amount of electric power and process 
steam that would be generated if the control values «(r)were applied to the 
actual power plant 3. The revenue term <?(r) represents revenues from sell- 

25 ing the electric power and process steam and is computed by multiplying the 
production with a price per unit of production. Computing the revenue q(r) 
preferably includes power prices that vary over time, either because agree- 
ments for providing a given amount of power at a given price and at a given 
time exist, or because estimates about future spot market power prices are 

30 available. 
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The costs c(r) are computed from the fuel useage, fuel costs, and process 
efficiencies required to produce, p{r) . The fuel consumption and process ef- 
ficiencies are computed using a plant simulation model and depend directly 
on the plant control values w(r) . Associating fuel costs with outputs is done, 
5 for example, by 

■ assigning all fuel costs to electric power output and considering the pro- 
duced steam as a free by-product, and 

■ assigning additional fuel costs to the steam output, if they arise because 
a higher steam production is required than would be produced at the re- 

io quired level of electrical power 

The rate of ageing e(r) represents costs associated with a reduction of life- 
time of plant components. Plant components are e.g. pressurized vessels, 
pipes or turbine blades. They are designed to withstand normal operating 

15 conditions for a given amount of time, i.e. their lifetime under nominal con- 
ditions. Subjecting a component to more severe operating conditions such as 
higher temperatures or temperature gradients increases component degra- 
dation and reduces component lifetime. This in turn increases long term 
costs for operating the power plant. According to the invention, an effect on 

20 these long term costs is included in the total operating costs. 

The rate of ageing e(r) is computed from models associated with various 
plant components, such as the following models, which are well know to one 
versed in the art: 

25 ■ Microcrack growth rate models for components operating at elevated 
temperatures, such as boiler pressure parts, steam pipes and steam tur- 
bine rotors. See for example "Life Assessment Technology for Fossil 
Power Plants", R. Viswanathan, Proc. American Power Conf., Vol. 57-111: 
1749-1768, 1995 and "A Lifetime Prediction Model for Thermomechanical 

30 Fatigue based on Microcrack Propagation", M.P.Miller et.al, ASTM STP 
1 1 86, 1993, pp. 35-49. 

■ Corrosion, oxidation, and erosion models based on extensions to mi- 
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crocrack growth models as described in "Corrosion Fatigue: Mechanics, 
Metallurgy, Elctrochemistry, and Engineering", T. Crooker and B. Leis edi- 
tors, ASTM 1983. 

■ Mechanical and thermomechanic stress models to determine stresses cre- 
ated in plant components, e.g. stresses in a steam turbine rotor shaft due 
to load oscillations and/or steam temperature changes. 

The rate of ageing e(r) represents the plant lifetime consumption or an asset 
value degradation expressed as a financial value under operating conditions 
caused by the control values u(j) . The entries in e(r) correspond to the ag- 
ing rate for each plant component and are, for example, computed directly 
from 

R da 

where a cnt is a predefined critical crack length for said component, a is the 
crack length computed using the models described above and dependent on 
w(r) , and, /?, is the component, cost, repair, or maintenance cost. 

In the inventive method, the following steps are performed 

a) The optimization module 1 determines the trajectory of plant control 
values w(r) from the present time / for a predetermined duration T and 
transmits said plant control values w(r)to the simulation module 2. The 
first time that this step is executed, plant control values are e.g. set to 
constant values for the duration T or are determined from stored stan- 
dard trajectories corresponding to a given operating condition such as a 
plant startup or a load change. Later executions of this step preferably 
use a variation of the plant control values w(r) of an earlier step. 

b) The plant model incorporated in the simulation module 2 determines 
from the trajectory of control values w(r), by simulation of the power 
plant behavior, a trajectory of the state of the power plant. From said 
state, corresponding trajectories of costs c(r) for operating the plant 
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and trajectories of revenues g(r) from selling produced output are de- 
termined, as well as the rate of ageing e(z) of plant components. In a 
preferred embodiment of the invention, trajectories of plant output val- 
ues p(r) are determined as well. 

A total operating cost J[u] is computed as an integral of an objective 
function that comprises the rate of ageing e(r) and a difference be- 
tween costs c(r) and revenues q{r). In a preferred embodiment of the 
invention, this total operating cost J[u] is defined as the functional 

t+T 

J[u]= J*,e(r) + * 3 [c(r)-$(r)] dr 



where R U R 3 are weighting matrices. The weighting matrices determine 
a relative influence of the magnitude of each component of each of the 
different vectors #(r),c(r),e(r). 

In a preferred embodiment of the invention, the above objective func- 
15 tion also includes a difference between plant output values p(r) and 

demand values d(r) . In this embodiment, the total operating cost J[u] 
is preferably defined as 

t+T 

J[u]= J *,e(r) + tf 2 p(r)-^r)f + * 3 [c(r)-*(r)] dr , 

t 

where R U R 2 ,R 3 are weighting matrices. The terms in the objective func- 
20 tion are interpreted as follows. The u d(r)-p(r) n term represents a pen- 

alty associated with not satisfying the required power output. The 
M c(r)-qr(r)" term represents costs minus revenue, i.e., an economic per- 
formance indicator, 
d) Steps a) through c) are repeated iteratively, with the optimization mod- 
25 ule varying the trajectory of plant control values u :[/,/ + T] i-> R n until an 

optimized lower total operating cost J[u] is arrived at. As the functions 
pW^(r),c(r),e(r) depend upon plant control values w(r), the optimiza- 
tion routine finds a strategy u :[t,t + 7]h R n that minimizes the total 
operating cost J[u]. 
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Examples for units are, 

■ for the rate of ageing e : $ per second 

■ for costs c or revenue q : $ per kW or $ per kg steam per second 
5 ■ for Rl : dimensionless 

■ for R2: $/kWA2 or $(s/kg)A2. 

■ for R3: kW or kg per second 

■ for demand d or production p : kW or kg per second 

10 In another preferred embodiment of the invention, the objective function is 

J[u]= J J?,e(r) + i? 2 Kr)-^)j+¥W-%"« dz , 

t 

where costs and revenues are weighted separately, using an additional 
weighting matrix R 4 , which allows for the cost and revenue vectors c(r), 
q(r) to have different dimensions. 

15 

The optimization step c) involves the problem of how to find the optimal 
solution efficiently, since the search space determined by the plurality of 
plant control values and the duration T is large. Standard optimization pro- 
cedures may be used, or dedicated solutions as described, for example, in 
20 US 5,347,466, US 5,886,895 or in "The Generalized Unit Commitment Prob- 
lem", Ross Baldick, IEEE Transactions on Power Systems, 1 0(1 ):465-475, Feb- 
ruary 1 995. 

In a preferred embodiment of the invention, at time /, the control strategy, 
25 i.e. the future values u :[t,t + T]\-+ R n corresponding to optimal plant opera- 
tion, is displayed graphically to an operator, providing him with guidance for 
controlling the plant. 

Figure 2 shows a block diagram corresponding to another preferred em- 
30 bodiment of the invention, in which, at time /, the control value u*(t) is ap- 
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plied to the actual power plant 3, resulting in automatic optimized plant op- 
eration. Preferably, measured plant values y(t) of the plant are input to the 
simulation module 2 in order to update the simulation model. Appropriate 
modeling and model update techniques are well known from the field opti- 
5 mal control . In particular, receding horizon control techniques might result 
useful. The power plant is, for example, a gas turbine, a combined cycle 
power plant, or of a boiler-turbine-generator type. 

A computer program product according to the invention comprises a com- 
10 puter readable medium, having thereon: computer program code means to 
make, when said program is loaded in a computer that is operationally con- 
nected to a power plant, the computer execute the method according to one 
of claims 1 to 6. A computer program according to the invention comprises 
said computer program code. 

15 

A power plant optimizing system that optimizes plant control values ac- 
cording to the invention comprises 

a) an optimization module 1 that is configured to determine a plurality of 
trajectories of plant control values w(r) and to determine a total operat- 

20 ing cost J[u] as an integral of an objective function that comprises a dif- 

ference between costs c(r) and revenues q(j) that are associated with 
operating the power plant according to the plant control values w(r) , and 

b) a simulation module 2 that is configured to determine, from the trajec- 
tory of plant control values w(r), by simulation of the power plant's be- 

25 havior, corresponding trajectories of costs c(r) for operating the plant 

and trajectories of revenues q(r) from selling produced output. 
The simulation module 2 is configured to determine, from the trajectory of 
plant control values u(r) , a trajectory representing a rate of ageing e(r) of 
plant components. The objective function comprises said rate of ageing e(r). 



30 



In a preferred embodiment of the power plant optimizing system, the simu- 
lation module 2 is configured to determine, from the trajectory of plant con- 
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trol values w(r), plant output values p{r), and the objective function com- 
prises a difference between plant output values p(j) and associated demand 
values d(r). 

5 The invention allows the operator or an automatic controller to incorporate 
not only short term financial considerations into a plant operation strategy, 
but also long term financial considerations that involve maintenance and re- 
placement costs and their dependency on operating conditions. 
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PATENT CLAIMS 



1. A method for optimizing power plant control values, where the method 
comprises the steps of 

5 a) determining a trajectory of plant control values w(r) , 

b) determining from the trajectory of plant control values i7(r) , by 
simulation of the power plant's behavior, corresponding trajectories 
of costs c(r) for operating the plant and trajectories of revenues 
g(r) from selling produced output, 

10 c) computing a total operating cost J[u] as an integral of an objective 

function that comprises a difference between costs c(r) and reve- 
nues q(r), 

d) iteratively repeating steps a) through c) with the optimization mod- 
ule varying the trajectory of plant control values u(r) until an opti- 
15 mized lower total operating cost J[u] is arrived at, 

characterized in that in step b), by modeling of plant components, a tra- 
jectory representing a rate of ageing e(r) of plant components is deter- 
mined from the trajectory of plant control values i7(r), and that, in step 

c) , the objective function comprises said rate of ageing e(r). 

20 

2. Method according to claim 1, characterized in that, in step b) trajectories 
of plant output values p{r) are determined, and that, in step c) the objec- 
tive function comprises a difference between plant output values p(r) 
and associated demand values d(r) . 

25 

3. Method according to claim 1, characterized in that the total operating 
cost J[u] is computed as 

/+r 

J[w]= J^e(r) + J? 3 [c(r)-9(r)] dr , 

where R U R 3 are weighting matrices, / is a present time and T is a prede- 
30 termined duration. 
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4. Method according to claim 2, characterized in that the total operating 
cost J[u] is computed as 

t+T 

J[u] = J R t e(r) + R 2 p(r) - p(r)f + * 3 [c(r) - *(r)] rfr , 

5 where R U R 2 ,R 3 are weighting matrices, t is a present time and T is a 

predetermined duration. 

5. Method according to claim 1, characterized in that at least part of the 
trajectory of plant control values «(r) corresponding to the optimized 

10 lower total operating cost J[u] is displayed to an operator. 

6. Method according to claim 1, characterized in that at least part of the 
trajectory of plant control values w(r) corresponding to the optimized 
lower total operating cost J[u] is used to control the actual plant 3. 

15 

7. A computer program product comprising a computer readable medium, 
having thereon: computer program code means to make, when said pro- 
gram is loaded in a computer that is operationally connected to a power 
plant, the computer execute the method according to one of claims 1 to 

20 6. 

8. A power plant optimizing system that optimizes plant control values and 
comprises 

a) an optimization module (1) that is configured to determine a plural- 
25 ity of trajectories of plant control values u(r) and to determine a to- 

tal operating cost J[u] as an integral of an objective function that 
comprises a difference between costs c(r) and revenues q(r) that 
are associated with operating the power plant according to the plant 
control values w(r), 

30 b) a simulation module (2) that is configured to determine, from the 
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trajectory of plant control values w(r), by simulation of the power 
plant's behavior, corresponding trajectories of costs c(r) for oper- 
ating the plant and trajectories of revenues g(r) from selling pro- 
duced output, 

5 characterized in that the simulation module (2) is configured to deter- 

mine, from the trajectory of plant control values w(r), a trajectory repre- 
senting a rate of ageing e(r) of plant components, and that the objective 
function comprises said rate of ageing e(r). 

10 9. A power plant optimizing system according to claim 8, characterized in 
that the simulation module (2) is configured to determine, from the tra- 
jectory of plant control values w(r), plant output values p(r), and that 
the objective function comprises a difference between plant output val- 
ues p(r) and associated demand values d{r) . 

15 
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ABSTRACT 

In a method and computer program product for optimizing power plant con- 
trol values and a power plant optimizing system an optimization module (1) 
minimizes total plant operation costs while achieving predetermined required 
5 output values for produced power and process steam. This is done by, at a 
given time, determining future values of control values and simulating, in a 
simulation module (2), the behavior of the plant up to a given future time. 
Corresponding fuel costs and generated power are determined in the simula- 
tion, and incorporated in an objective function. The optimization module (1) 
10 minimizes the objective function by varying the control values. According to 
the invention, a rate of ageing of plant components is determined when 
simulating the future behavior of the plant, and the objective function to be 
minimized comprises said rate of ageing. 

15 (figure 1) 
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